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We show how an analysis of the multiple muon events in 
the Monopole Astrophysics and Cosmic Ray Observatory detector 
can Be used to determine the cosmic ray composition. 
Particular emphasis is placed on the interesting primary 
cosmic ray energy region above 2000 TeV/nucleus. 

An extensive study of muon production in cosmic ray 
showers has been done by Gaisser and Stanov (2). Results 
were used to parameterize the characteristics of muon 
penetration into the earth to the location of a detector. 

The mean number of muons at et slant depth X is 

<N > - | (.0142)(^£)- 775 (1 - ^) 5 ' 96 sec0, 

^ fJ V P 

where E is the primary energy per nucleon, 9 is the primary 
direction relative to the vertical, and 


E^(TeV) - . 53 ( e ’ 4x -l ) 


is an effective muon threshold energy (so that E > E is 
required). Each nucleon in a primary nucleus contributes 
this same average number of muons. The N distribution is 
taken to be Poisson. ^ 

The lateral distribution of muons at the detector depth 
is 


dN 


dR 


— a R e , where 


<v 


3.13E * 46 + 13. 2E • 31 (=^)' 62 

U u 'E ' 


Here R is measured perpendicular to the shower core. The R 
and Maldistributions are uncorrelated. ^ 

These parameterizations yield predictions for the Kolar 
Gold Field experiment (3) which agree well with the data (4). 
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Fig. 1 a) Distribution of overburden 
(depth) that muons must penetrate to 
reach the detector, b) Distribution 
of primary cosmic ray energies (per 
nucleus) for all depths and for 
depths greater than 7 KMWE. 

For example, in Fig. lb we show the d 
energies (per nucleus) for all depths 
depths that are greater than 7 KMWE. 
events with primaries < 2000 TeV/nucl 
events but only 8% for the restricted 


The MACRO detector 
is taken to be a 
horizontal rectangle 
12m by 112m, and all 
muons striking this 
rectangle are assumed 
to be detected. (We 
are currently 
redesigning the 
apparatus so that it 
will have a significant 
size in the vertical 
direction; the present 
height is 5.7 m high) . 

Muons that reach 
the detector pass 
through a variable 
slant depth of material 
in the 

irregularly-shaped Gran 
Sasso Mountain, in 
Fig. la we show the 
slant depth 
distribution in units 
of KMWE. Because we 
measure muon directions 
with Q high resolution 
(~.2 ), we can select 
on the slant depth and 
thereby have a data 
sample with a 
distribution of primary 
energies that can be 
changed. 

istribution of primary 
, and separately for 
The percentage of 
eus is 31% for all 
sample . 


A major goal of MACRO is to determine the cosmic rav 
composition above the 2000 TeV "knee" in the energy Y 
distributions. We have studied four models for this hiqh 
energy composition: y 


1. Fe All primaries are iron nuclei. 

2. Md A mix of primaries is used with iron dominating 

at high energy. 3 

3. LE Primary composition at high energies is the same 

as at 50 GeV/nucleon. 

4. p All primaries are protons. 

Models Md and LE are described in reference 3. 
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In Fig. 2 we show the detected muon multiplicity 
distributions for the four models (solid curves), all 
normalized to the single muon rate. 

For a given nucleus energy, the proton has a much larger 
energy per nucleon than iron, so protons give many more 
single muon events; however, the iron nucleus produces many 
more muons so its multiple muon to single muon ratio is 
larger. However, lateral spreads (decoherence) from iron are 
larger, so a large detector like MACRO is needed to exploit 
this difference. 



MUON MULTIPLICITY, NM 


Fig. 2 Muon multiplicity distri- 
bution divided by single muon rate 
for various hypotheses. (See text.) 


The dashed curves in Fig. 2 show how difficult it is for 
a small detector (10m* area) to measure high multiplicities. 
For 10 fj’ s, the Fe hypothesis yield is down two orders of 
magnitude. The absolute MACRO yields are about 500K/year for 
Fe and 13M/year for protons, including single muon events. 
Depending on composition, the number of Events with NM > 5 
will be 500 to 20,000! 
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One aspect of the muon decoherence which is particularly 
sensitive to primary composition is the maximum separation of 
a pair of muons in a multiple muon event. In Fig. 3 we show 
this distribution for the four composition hypotheses and for 
three different ranges of slant depth (and thus three 
different primary energy distributions). 

Figures 2 and 3 both show that MACRO is sensitive to the 
cosmic ray primary composition at high energy. 



Fig, 3 Distribution of maximum two- 
muon separation in a multiple muon 
event for three ranges of overburden. 
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